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Abstract 

The spectra of lepton pairs from correlated open charm and bottom decays in ultrarel- 
ativistic heavy-ion colhsions are calculated. Our approach includes energy loss effects of 
heavy quarks in deconfined matter which are determined by temperature and density of 
the expanding parton medium. We find a noticeable suppression of the initial transverse 
momentum spectrum of heavy quarks due to the energy loss at LHC conditions. Within 
the central rapidity covered by the ALICE detector system, the dominant contribution 
from bottom decays to the high invariant mass dielectron spectrum is predicted. 
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The production of high invariant mass dileptons in relativistic heavy-ion colhsions has 
attracted a great deal of interest in the recent years from both the experimental and 
theoretical sides. Dilepton measurements are planned with the PHENIX and ALICE de- 
tector facilities at the relativistic heavy-ion collider RHIC in Brookhaven and at the large 
hadron collider LHC in CERN, respectively. In despite of many theoretical attempts the 
nature of the dominating source of dileptons, which will be measured in these detectors, 
is still unclear. Among the candidates that can give strong contributions in the large 
invariant mass continuum region near and beyond the J / ip one can consider the following 
processes in which dileptons are generated: (i) hard initial quark - anti-quark collisions 
give rise to the Drell-Yan yield, (ii) semi-leptonic decays of heavy quarks like the charm 
and bottom, which are also produced dominantly in hard initial parton interactions, and 
(iii) the thermal dilepton radiation resulting in interactions between secondary partons or 
even in a locally thermalized quark-gluon plasma (QGP). The yield from (iii) is expected 
to serve as a direct probe from the QGP [|T], 0]. 

It is commonly believed that, due to the large invariant mass M, at least the Drell-Yan 
process and the initial heavy quark production in pp and pA collisions are under reliable 
theoretical control by means of perturbative QCD ^, |] and, therefore, can provide a 
reference for the thermal dilepton signal from deconfined matter. At the same time all up 
to date calculations of high-M thermal dilepton rates from early parton matter contain 
uncertainties due to the lack of knowledge of precise initial condition parameters such as 
the energy density and quark-gluon phase space saturation. The parton cascade model 
and calculations based on the HIJING event generator differ substantially in predicting 
the thermal high-M dilepton spectra. Recent estimates within the mini-jet mechanism 
of quark-gluon matter formation p[, with initial conditions similar to the self-screened 



parton cascade model |T0[ , point to a strong competition between the Drell-Yan yield and 
the thermal signal in a wide range of beam energies of the colliding nuclei. 

With respect to the wanted dilepton signal from the QGP the dileptons from charm 
and bottom decays play a role of substantial background The magnitude and the 
characteristic properties of such background is presently matter of debate. As shown 
by Shuryak [|12| for RHIC conditions the energy loss effects of heavy quarks propagating 
through a medium can drastically reduce the dilepton yield from charm and bottom decays 
in the region of large invariant dilepton masses. Further calculations [jl3[, including the 
PHENIX acceptance at RHIC, also predict a very strong suppression of dileptons from 
heavy quark decays, which then is partially below the Drell-Yan yield. The studies in 
refs. [Il2| , [T3[| are based on the approximation of a constant rate of energy loss —dE/dx = 1 



or 2 GeV/fm. Independently of the assumed mechanism of energy loss the approximation 
of a constant value of dE/dx is difficult to justify for strongly expanding and cooling 
matter. In general, the value of dE/dx is expected to depend not only on the geometry 
and a fixed mean free path of heavy quarks but also on the temporal evolution of the 
medium. As shown below, in particular within the gluon radiation mechanism of energy 
loss Ijl^, |15|, |T^, dE/dx can depend on widely varying characteristic quantities of the 



parton medium formed at RHIC and LHC 

In the present paper we study the dileptons resulting from decays of correlated open 
charm and bottom produced initially in hard collisions and undergoing an energy loss 
which is determined by the temperature and the density evolution of the deconfined 



medium. Going beyond refs. [|T2|, |T3[ we concentrate here on LHC conditions, where the 
energy loss appears to be stronger as compared to RHIC. In a certain scenario this leads 
to almost complete stopping of initial charm and bottom in deconfined matter produced 
at LHC. The central rapidity acceptance of ALICE will select in this case the electron 
pairs mainly from bottom decays. It is assumed that, due to a like-sign subtraction, 
the combinatoric background from uncorrelated pairs can be removed, and henceforth we 
focus on correlated pairs. 

In the framework of perturbative QCD the number of the heavy quark - anti-quark 
(QQ) pairs, produced initially with transverse momenta p±i = —p±2 = Pi_ at rapidities 
yi and y2, in central AA collisions can be calculated from 

dNqQ = T4a(0) ICq H{yi, y2,p±) dp]_ dyi dy2, (1) 

where H{yi,y2,p±) is the standard combination of structure functions and elementary 
cross sections (see for details |11|) 

f da^ 

H{yi,y2,p±) = xiX2<fg{xi,Q'^)fg{x2,Q'^)-^ 

+ J2[f,{x„Q')Mx2,Q') + f,{x2,Q')fg{x^,Q')]^^, (2) 
where /i(a;, Q^) with i = g,q,(j are the parton structure functions, Xi^2 = 



m_|_ (exp{±yi} + exp{±y2}) / \fs and mj_ = yp^. + "^q- heavy quark masses we take 
nic = 1.5 GeV and = 4.5 GeV. We employ throughout the present paper the HERA 
supported structure function set MRS D'- |jl8| from the PDFLIB at CERN. The overlap 
function for central colhsions is TaaIO) = A'^/{TrR\) with Ra = I.IA^/^ fm and A = 200 
in this paper. We here do not include shadowing effects of nuclear parton distributions 
since for heavy quark production they are expected to be not very important and can 



be considered separately. Our calculation procedure is based on the lowest-order QCD 
cross sections da^gjdt for the subprocesses gg QQ and qq QQ with the simulation 
of higher order corrections by the corresponding K,q factor. Such a procedure reproduces 
within the needed accuracy the next-to-leading order calculations of the heavy quark 
pair distributions with respect to their invariant mass, total pair rapidity and rapidity 
gap. We find the scale = 'imq and Kq = 2 as most appropriate. 

For the Drell-Yan production process of leptons at rapidities yi and j/2 and transverse 
momenta p^i = —p±_2 = P± we have 

dN^^ = TAAiO)JCDYL{yi, 1/2, p±)dp]_dyidy2, (3) 

r ^ ^ 1 da'-'' 

L{yuy2,p±) = Xl^i^^'^^'^^i'^^) + /9(a^2,Q^) /g-(a;i,(5^)J -4-, (4) 

q,q "''^ 



with da'^/di = ^ch(yi - ^2)/ (2pl[l + ch(?/i - 2/2)]^), 2:1,2 = P± (exp{±?/i} + exp{±?/2}) 
/v^, a = 1/137, = X1X2S = M2 and JCdy = 1-1 1- 

Since the energy loss of heavy quarks depends on the properties of produced matter, 
which is propagated through, one has to specify the initial conditions and the space-time 
evolution of such a medium. To do this in a coherent manner we employ for the mini- 
jet production (which dominates the parton matter formation) the same lowest-order 
approximation and parton structure functions as for the initial heavy quark production. 
Adding a suitably parametrized soft component 0] we get the initial temperature Tj = 
1000 (550) MeV, gluon fugacity Af = 0.5 and light quark fugacity = Af/5 of the mini- 
jet plasma formed at initial time Tj = 0.2 fm/c at LHC (RHIC). The space-time evolution 
of the produced parton matter after is governed by the longitudinal scaling-invariant 



expansion accompanied by quark and gluon chemical equilibration processes |19|, ^ . We 



take for definiteness full saturation, i.e. A^ = A^ = 1 at deconfinement temperature 
Tc = 170 MeV. (Due to gluon multiplication this seems to be justified, while the 
chemical equilibration of the quark component might be somewhat overestimated by this 
ansatz.) The actual time evolution is determined by e oc r~^/^ for energy density and 
\q,g Q(- -^2 ipg^ which must be inverted to get T(r) explicitly. This results in final pion 
rapidity densities (see [^) which are in agreement with other estimates 0. 

To model the energy loss effects of heavy quarks in expanding matter we assume 
as usual in the Bjorken scaling picture that a heavy quark produced initially at given 
rapidity will follow the longitudinal flow with the same rapidity. Therefore, with respect 
to the fluid's local rest frame the heavy quark has essentially only a transverse momentum 
p± which may depend on the proper local time r in accordance with the energy loss in 
transverse direction. (Note that the transverse expansion at early times can be neglected. 



and we take as averaged transverse radius of parton matter R = 7 fm.) 

Recently the total energy loss of a high-energy parton propagating transversally 
through expanding deconfined matter has been derived as (i-E/(ia;| expanding = C,dE/dx\Tf 
ITf , supposed the medium cools according to a power law (or similar) and the initial 



time tends to zero. The numerical factor ^ is of the order 2 (6) for a parton created 
inside (outside) the medium. The subscript Tj indicates that the stopping power depends 
only on the final temperature at which the deconfined medium is left. We shall denote 
hereafter this stopping scenario as model I and take dE/dx\Tj from |l^, |l^, [l^ with 
the new, improved numerical factors. 

We are going to contrast model I with another one which relies on an accumulative 



energy loss similar to [0, |T3|. To calculate the evolution of the transverse momentum 
p_i_(r) of quarks propagating the distance ^^(r) in transverse direction we adopt the results 
of fl^ for the energy loss of a fast quark in a hot QCD medium (model II) 



dpi. Aash 2 , 2N1/4 1 ( yP± + ^Q 1 /^x 

^ = -371^^" + "'^^ ^"1, Lfcg I' 
dr± 



dr + ml 



(6) 

where the parameter = 2mlf^{T) = |(Ag + ^Xq)nasT^ fl^, ^ is used as an average 
of the momentum transferred in heavy quark-parton scatterings, and the strong coupling 
is described by as = 0.3. For the mean free path of gluons we take L ^ {o'ggfT'g)^^ 
with integrated cross section agg ~ a^/kl and gluon density Ug ~ AgT^, resulting in 
= 2.2 OgT [|T9| , Strictly speaking, eqs. (5,6) are valid for a static medium under 



certain conditions [|T5|. Integration of these equations therefore neglects the destructive 
interferences of the Landau-Pomeranchuk-Migdal effect. The r.h.s. of eq. (^ is actually 
the rate of energy loss dE/dx which is essentially oc (p^ + rnQY^^. Due to the expansion 
of the matter dE/dx depends on the time. For instance, at LHC conditions and initial 
momentum of a heavy quark of 10 GeV the stopping power derived from eq. (5) drops 
from 6 GeV/fm down to 0.3 GeV/fm during the expansion of deconfined matter. 

To get the spectra of charm and bottom quarks after energy loss we use a Monte 
Carlo simulation with a uniform distribution of the random initial position and random 
orientation of QQ pairs in the transverse plane. In model I the total energy loss is 
determined by the transverse distance ci of a created heavy quark to the boundary of the 
system, the quark initial energy and the temperature when leaving the system. Most 
quarks experience an energy loss oc d and only a few ones oc d"^. 

In model II we integrate eqs. (||,|) together with the above described time evolution 



of T(r) and Ag,c,(r). Eq. is used to check whether the considered heavy quark propa- 
gates still within deconfined matter; if it leaves the deconfined medium it does not longer 
experience the energy loss according to eq. (^). (The subsequent energy loss in a mixed 
and hadron phase might be incorporated along the approach of [l^-) A considerable 
part of heavy quarks can not escape the parton system thereby undergoing thermaliza- 
tion. Following we consider a heavy quark to be thermalized if its transverse mass 
m_L during energy loss becomes less than the averaged transverse mass of thermalized 
light partons at given temperature. These heavy quarks we assume to be distributed as 
dNQ/dm± oc m\exp{—m±/Tcs} with an effective temperature Tefj = 150 MeV. 

The results of our calculation for transverse momentum spectra of charm and bottom 
quarks at midrapidity are shown in fig. 1. Within model II we obtain a strong suppression 
of the initial transverse momentum spectrum at large transverse momenta due to the 
energy loss both for charm and bottom. The bump below 3 GeV is caused by the thermal 
redistribution. In comparison with RHIC conditions [|T^, ^ the suppression at LHC is 
found to be more pronounced due to the longer life time of deconfined matter and its 
higher initial temperature and density. It should be stressed that such a drastic change of 
the p± spectrum of heavy quarks in AA collisions, in comparison with pp or pA collisions, 
would be a measurable effect which for itself can help to identify the creation of a hot and 
dense parton gas ^ Q . In contrast to model II, model I results in a smaller depletion of 
the large-p^ spectrum (see fig. 1). 

Basing on the analysis of ref. |]ll| we employ a delta function like fragmentation func- 
tion for heavy quark fragmentation into D and B mesons. It gives the same transverse 
momentum for the meson as for the parent heavy quark. The correlated lepton pairs from 
DD and BB decays have been obtained from a suitable Monte Carlo code which employs 
the single electron distributions as delivered by JETSET 7.4. The average branching ratio 
of -D — i> / + X is taken to be 12%. We consider here only the so-called primary leptons 
[|13| which are directly produced in the bottom decay b I X with the average branching 
ratio 10%. The secondary leptons I' produced in further branchings like b —>■ cX —>■ I'Y 
have much smaller energy and do not contribute to the large invariant mass region. Such 
simplification does not affect significantly the dielectron yield at M > 2 GeV as we have 
checked by comparison with ref. [jlll. To get the dielectron spectra from charm and bot- 
tom decays at LHC we take into account the acceptance of the ALICE detector: the 



electron pseudo-rapidity is \rje\ < 0.9 |25] with an additional cut on the electron trans- 
verse momentum p\> 1 GeV to reduce the lepton background from light hadron decays. 
We do not take into account the limited high-p_|_ acceptance (p^ < 2.5 GeV [^) in our 



simulations to demonstrate the general tendency of the spectrum in the high-M region. 
Our results of heavy quark decays at LHC are depicted in fig. 2 for charm (a) and bottom 
(b) for M > 2 GeV. Due to the strong energy loss in model II the contribution from charm 
decays drops below the Drell-Yan yield. At the same time the contribution from bottom 
decays becomes dominant up to M = 6 GeV. This can be understood also on a qualitative 
level. Even after strong energy loss the massive B mesons have enough energy to produce 
electrons with p']_ > 1 GeV. Additionally the central rapidity cut of the ALICE detector 
selects electrons from charm being near mid-rapidity, while for bottom the corresponding 
rapidity region is greater. Further, we actually estimate a lower limit for bottom decay; 
as discussed above the secondary electrons contribute somewhat at not too large values of 
M. Therefore, at LHC the dielectron continuum region above 2 up to 6 GeV is expected 
to be dominated by bottom decays. This prediction is also supported by model I. Here 
the energy loss is smaller, and the lepton pairs from bottom dominate at larger invariant 
mass. At M < 5 GeV bottom and charm are equally strong, however an additional energy 
shift, e.g. caused by including a non-5 like fragmentation, lets bottom win. (The Peterson 
fragmentation function results in a > 3 times stronger bottom contribution in comparison 
with charm. Also a slightly changed numerical factor, say ^ = 3, suppresses the charm 
contribution down to the Drell-Yan yield, while the bottom contribution is less affected). 

To demonstrate the general tendency of the relative contribution from the Drell-Yan 
process and heavy quark decays with the change of the collider energy we plot in fig. 3 
the results of our calculations for RHIC initial conditions including the PHENIX detector 
acceptance for the electrons, i.e. \rie\ < 0.35 |2^. As above for LHC we here also employ 
the gate > 1 GeV. For model II one observes in fig. 3 a strong competition of the 
Drell-Yan process and correlated heavy quark decays in a wide region above M = 2 GeV 
similar to |T^. Imposing an additional kinematical cut, which selects back-to-back pairs, 
Drell-Yan dileptons very clearly dominate. However, notice that model I does not give a 
noticeable effect of energy loss at RHIC conditions. 

In summary, taking into account the energy loss effects of heavy quarks we calculate 
the spectra of electron pairs from correlated open charm and bottom decays at LHC con- 
ditions within two different scenarios. Within the naive scenario with accumulative energy 
loss in the expanding deconfined matter we obtain a measurable, strong suppression of 
large transverse momenta of heavy quarks, while the recent treatment of the energy loss 
1^ predicts a smaller effect. With respect to these different results and uncertainties 



m 



inherent in our approach (like the neglect of quark masses in some energy loss formu- 
lae, finite energies of quarks, use of leading order results in the strong coupling constant. 



neglect of the confinement transition and the influence of confined matter) any predic- 
tion seems at best semi-quantitative. Therefore, an experimental determination of the 
transverse momentum spectra of heavy quarks in future heavy-ion collisions appears as 
challenging task to pin down the essential effects. Qualitatively, our calculations predict 
the dominating contribution of bottom decays in the high invariant mass region of electron 
pairs which are planned to be measured at central rapidity in the ALICE detector. 
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Figure captions 



Fig. 1: The p± spectrum of heavy quarks ((a): charm, (b): bottom) at midrapidity in 
AA coUisions at LHC (a/s = 5500 AGeV). The dotted curves depict the initial production 
without energy loss, while the dashed and solid lines show our result with energy loss 
according to model I and II and Teflf = 150 MeV for thermahzed heavy quarks. 

Fig. 2: The invariant mass spectrum of dielectrons from correlated charm (a) and 
bottom (b) decays as well as Drell-Yan electron pairs filtered throughout the LHC-ALICE 
acceptance, but the high-pj_ cut of 2.5 GeV is here not implemented. 

Fig. 3: The invariant mass spectrum of dielectrons from correlated charm (dotted line) 
and bottom (sohd line) decays as well as Drell-Yan electron pairs (dashed line) filtered 
throughout the PHENIX acceptance. The energy loss is calculated within model II. 
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